High-energetic Cosmic Antiprotons from 
Kaluza-Klein Dark Matter 



Torsten Bringmann 

Department of Physics, AlbaNova University Center, Stockholm University, SE - 106 
91 Stockholm, Sweden 

E-mail: troms@physto.se 

Abstract. The lightest Kaluza-Klein particle (LKP) in models with universal extra 
dimensions is an interesting dark matter candidate that has recently received great 
attention. Here, we investigate the antiproton flux from LKP annihilations in the 
galactic halo. In our analysis we include different halo density profiles and allow for 
part of the dark matter to be concentrated in 'clumps' rather than being distributed 
homogeneously. After re-analyzing the observational bounds on the allowed amount 
of dumpiness, we find that LKP annihilations may well give a significant contribution 
to the antiproton flux at energies higher than about 10 GeV, while for energies above 
around 500 GeV the conventional background is expected to dominate again. The 
shortly upcoming PAMELA satellite will already be able to measure part of this high- 
energy window, while planned experiments like AMS-02 will have access to the full 
energy range of interest. 
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1. Introduction 

The observational evidence for the existence of dark matter has become overwhelming 
and covers an impressive range of distance scales, from galaxy rotation curves [T] via the 
mass-to-hght ratio of galaxy clusters j2] to large scale structure surveys jH] . Adding our 
knowledge from microwave background observations |1] , supernova luminosity distances 
[S] and big bang nucleosynthesis (BBN), one can determine the dark matter density of 
the universe to a remarkable accuracy. Yet, the nature of the dark matter still remains 
a mystery. Weakly interacting massive particles (WIMPs), however, are particularly 
attractive dark matter candidates, as gauge couplings and masses at the electroweak 
symmetry breaking scale give a thermal relic abundance that naturally comes within 
the observed range p. 

A well-studied example of such a WIMP candidate is the neutralino, the lightest 
supersymmetric particle. We concentrate instead on an interesting alternative that 
has received much attention in recent years - Kaluza-Klein (KK) dark matter that 
appears in models with large extra dimensions [313 EI- In contrast to the neutralino, 
which is a Majorana fermion, this dark matter candidate is a massive vector particle, 
and consequently there is a rich, new phenomenology to explore that may well apply 
also to other dark matter candidates that are non-Majorana in nature. Observational 
prospects have been explored in some detail ^Ij and found to be promising, both for 
direct-detection experiments with very sensitive detectors ^T], as well as for indirect 
detection through a contribution to the cosmic ray flux in photons [121 Uni, positrons 
[Ti] and neutrinos ^3] of all energies up to the dark matter particle's mass. 

Here, we present the first analysis of the expected contribution to the antiproton 
flux. Compared to supersymmetry fOl 1^ , the a priori interesting feature of KK 
dark matter in this respect is a rather large branching ratio into quarks and gluons, and 
a mass at the TeV scale, which is higher than what is usually the case for neutralinos. 
Consequently, the focus of our attention will be on a possible distortion of the antiproton 
spectrum for high energies (above several GeV). Traditionally, one has instead often 
focused on low energies (below and around the peak at about 1 GeV) when looking for 
exotic physics in the antiproton spectrum, which has the disadvantage of being plagued 
by much greater uncertainties on the theoretical side. 

This paper is organized as follows. We begin by introducing the model that gives 
rise to the KK dark matter candidate B^^\ In Section IHl we then present the antiproton 
yield from B'^^'> annihilations in the galactic halo. Since the expected flux depends 
to a large extent on the details of the dark matter distribution, different halo profiles 
will be considered in Section jH This includes a detailed discussion of the possible 
existence of dark matter substructure ("clumps"), with updated bounds on the amount 
of dumpiness that is observationally allowed in these kinds of scenarios. The two-zone 
cylindrical diffusion model that is used to propagate the antiprotons through the galaxy 
is presented in Section O and in Section IHl the expected antiproton spectrum is shown 
and confronted with observational prospects. Finally, Section |7| concludes. 
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2. Universal extra dimensions 

The lightest KK particle (LKP) that appears in models of universal extra dimensions 
(UED) [7], is the first viable particle dark matter candidate to arise from extra- 
dimensional extensions of the standard model (SM) of particle physics (for the first 
proposal of TeV sized extra dimensions, see jlH])- In UED models, all SM fields 
propagate in the higher-dimensional bulk. For the effective four- dimensional theory, 
this means that every particle is accompanied by a tower of increasingly more massive 
KK states. Momentum conservation in the extra dimensions, which leads to KK mode 
number conservation in the effective four dimensional theory, is broken at the boundaries 
of the compact space by an orbifold compactification that projects out unwanted degrees 
of freedom at the zero mode level. If the boundary terms introduced at the orbifold fixed 
points are identical, a remnant of KK mode number n conservation is left in the form of 
conserved KK parity, (—1)", and the LKP is stable |H]. This is analogous to conserved 
R-parity in supersymmetric models, which in that case ensures the stability of the 
lightest supersymmetric particle. An attractive feature of this dark matter candidate, 
as compared to the supersymmetric neutralino, is that the unknown parameter space is 
quite small and will be scanned throughout by next generation's accelerator experiments. 

We will consider the simplest, five dimensional model with one UED compactified 
on an /Z2 orbifold of radius R. At tree-level, the nth KK mode mass is then given 



where msw is the zero mode mass of the corresponding SM particle. However, the 
identification of the LKP is nontrivial because radiative corrections to the mass spectrum 
of the first KK level are typically larger than the corresponding electroweak mass shifts. 
A one-loop calculation shows that the LKP is well approximated by the first KK mode 
of the weak hypercharge gauge boson B^^^ [8J. The i?*^^^ relic density was determined 
in jni EOj- Depending on the exact form of the mass spectrum and the resulting 
coannihilation channels, the measurements from the Wilkinson Microwave Anisotropy 
Probe (WMAP) |1| of ficDM/i^ = 0.12 ± 0.03 corresponds to 0.5 TeV < m^d) < 1 TeV 
if all dark matter is made of LKPs. Here, ^cdm is the ratio of cold dark matter to 
critical density and h is the Hubble constant in units of 100 kms~^Mpc~^. Collider 
measurements of electroweak observables give a current constraint of R^^ > 0.3 TeV 
[71 121], whereas LHC should probe compactification radii up to 1.5 TeV j23- 

3. Antiprotons from LKP annihilations 

The number of antiprotons p per unit time, energy and volume element that are 
produced by LKP annihilations at some position r in the galactic halo is given by 
the source function 



by 




(1) 
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Table 1. Branching ratios (in percent) for all annihilation channels at tree level that 
are relevant for antiproton production. The values shown above are calculated for LKP 
masses mg(i) — 0.8 [0.5, 1.0] TeV and several values for the mass shift between the 
LKP mass and the mass M of the other first level KK modes. 



where T is the antiproton kinetic energy, p^ci) the B^^^ halo density and its mass. 

The sum runs over all annihilation channels /, where B^^ and dN^ /dT are the branching 
ratios and fragmentation functions into antiprotons, respectively. The factor 1/2 appears 
due to the fact that B^^^ particles always annihilate in pairs; note that this factor has 
often not been taken into account when calculating antiproton (or other cosmic ray) 
yields. Finally, the B^^^ annihilation rate is given by [Oj 

(aan„t;)B(i) - 3 ■ 10-26(0.8 TeV/m5(i))2 cm^ s"^ . (3) 

A considerable source of astrophysical uncertainty lies in the unknown properties 
of the dark matter halo profile pcdm (which in the following will be assumed to be 
dominated by LKPs, i.e. Pcdm — Pbw)i this will be discussed in detail in the next 
section. In this section, we will instead concentrate on the antiproton differential 
spectrum per annihilation event 

/ 

that can be computed to a much better accuracy. At tree-level, with all other first level 
KK modes degenerate in mass, B^^'' pairs annihilate into quark pairs (35%), charged 
lepton pairs (59%), neutrinos (4%), charged (1%) and neutral (0.5%) gauge bosons and 
Higgs bosons (0.5%). Table^shows the exact branching ratios for all final states that are 
relevant for antiproton production. Clearly, they are very insensitive to the parameters 
of the model. 

The annihilation into gluon pairs is loop-suppressed, but since one can expect a 
rather large antiproton yield from gluon hadronization, it seems worthwhile to include 
this annihilation channel as well. Recently, the similar process B^^^B^^^ 77 has been 
analysed in detail ^3]. For a two-gluon final state, only the fermion box diagrams 
appearing in that analysis contribute, with of course no contribution from leptons. The 
only potential contribution from loops containing scalars comes from diagrams with an 
s-channel Higgs coupled to a quark loop (to which also the final gluons couple); these 
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Figure 1. The differential antiproton spectrum per B'^^^ pair annihilation event. From 
bottom to top, the plots correspond to an LKP mass of = 500, 800 and 1000 

GeV, respectively. 

amplitudes, however, are suppressed by a factor niq/m^d) due to the Yukawa couphng 
and can thus be neglected if. One can therefore easily derive an exact expression, up to 
C(mg/m^(i)), for the total cross section of B^^^B^^^ — > gg: Just take Eq. (16) in |Tn] 
and substitute the electric charge Q with 1 (this corresponds to setting g^j-j = 11/2); 
the color average is then performed by the simple recipe ttg^ — > (2/9) [221 • The result 
is a cross section for gluon pair production that is roughly 0.5% of the total annihilation 
cross section at tree level, i.e. about the same size as the (tree-level!) contributions from 
massive vector bosons. 

We use the tabulated fragmentation functions dN^ /dT of the DarkSUSY package 
|24j . which are based on a large number of events (10^ per annihilation channel and 
mass) as simulated with the Monte Carlo program Pythia 6.115 [22] • (While the 
parametrizations given in reproduce the spectrum to a very good accuracy for low 
and intermediate energies, they overpredict the highest-energy tail by a factor of up to 
3). The fragmentation function for Higgs final states is not explicitly given; it can easily 
be obtained by letting the Higgs decay in flight and then boosting the fragmentation 
functions for those decay products that are relevant for antiproton production. The 
resulting total differential antiproton spectrum g{T) is shown in Fig. ^for different B^^^ 
masses. Analyzing the individual terms in the sum (jH) separately, one finds that top and 

J Doing the full analytical calculation it is found that, even for the top quark, this diagram contributes 
only 1.7% to the cross-section of the full process B^^^B'-^^ gg. 
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charm quarks give as expected the largest contribution, except for the highest energies, 
where the W^W~ channel starts to dominate. 

4. Halo profiles 

As described in the last section, the ingredients for the antiproton source function Qp 
that derive from the underlying particle physics can be computed in a straight-forward 
way to a satisfying degree of accuracy. Unfortunately, the same does not hold for the 
astrophysical part; the distribution of dark matter in our galaxy is to a large extent 
unknown and since pcDM enters Qp quadratically, this puts considerable hmits on the 
possibility of making reliable predictions for the expected flux of cosmic rays. 

4.1. Smooth dark matter distributions 

Direct observations can only poorly constrain the form of the halo profile in the Milky 
Way, so it is usually taken from N-body simulations of gravitational clustering. Since 
these simulations only reach resolutions of about 0.1 kpc for galaxies with the size of 
the Milky Way, the innermost slope of the density profile has to be extrapolated and is 
thus bound to a considerable amount of theoretical uncertainty. Having this in mind, a 
generic parameterization of possible spherical halo profiles is given by 



where the indices (a, /3, 7) fix the halo model in question, Rq = 8.5 kpc is the distance 
of the sun to the galactic center and po ~ 0.3 GeV cm~^ the local halo density, which 
is known only up to a factor of almost 2 j2n|. The scale radius a determines where 
the transition for the radial dependence of Pcdm from large (pcdm oc r~^) to small 
(pcDM oc r~'^) galactocentric distances takes place. N-body simulations favour cuspy 
halo distributions like the NFW profile |2ZI with (a,/?, 7) = (1,3, 1), or the even more 
cuspy Moore profile ^] with (1.5,3,1.5). In fact, there is also analytic support for 
considering profiles with 7 > 1 j^H] ■ Here, we will for comparison also take into account 
the rather conservative case of an isothermal sphere, (2, 2, 0), that has a constant density 
core. There are indications that the scale radius a is not an independent parameter but 
strongly correlated with the virial mass of the galaxy SHj. In the case of the Milky Way 
and the isothermal sphere (NFW, Moore) profile, one finds a = 4 kpc (21.7 kpc, 34.5 
kpc) [HI], which we will use in the following as reference values. 

The most recent simulations actually suggest a universal halo profile of the form 



with n ~ 5 |32j . However, taking into account the effects of baryons in a simple model 

of adiabatic contraction jHS], this profile is effectively transformed to a Moore profile 
[HI]. Therefore, we restrict ourselves in the following to the isothermal sphere, NFW 
and Moore profiles. This should cover a plausible range of realistic halo models, serving 




(5) 



PcDM(r) ^ Po exp - (r/ a) 



(6) 
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to illustrate the uncertainties that derive from our lack of knowledge of the exact form 
of the galactic dark matter distribution. 

4-2. Dark matter substructures 

In cosmologies with cold dark matter, structure forms hierarchically, so it is reasonable to 
assume the existence of substructures ( "clumps" ) within the smooth halo distributions 
discussed in the last section. This is very interesting from the point of view of cosmic 
rays as an indirect detection method of dark matter, since the total expected flux is 
always obtained by averaging Pcdm o'^^ way or the other - which generically gives 
larger values for inhomogeneous dark matter distributions. This effect has mainly been 
studied for substructures larger than lO^M©, roughly correponding to the resolution 
limit of the N-body simulations mentioned earlier. Recently however, attention has 
turned to much smaller scales when it was discovered that for every WIMP dark matter 
candidate there is a natural cut-off in the power spectrum, below which density 
perturbations are washed out by collisional damping and subsequent free streaming [32] • 
A considerable part of the dark matter should therefore have collapsed into substructures 
with a mass of about M^, thereby forming the first gravitationally bound objects in the 
cosmological evolution. Based on very recent, high-resolution N-body simulations, it was 
estimated that a relatively large fraction of these microhalos can survive until today, 
giving a contribution of several percent to the total amount of dark matter in the solar 
neighborhood [SE] • These results are based on the assumption of neutralino dark matter, 
so eventually the whole analysis has to be redone for the case of KK dark matter |37j ; 
while this can certainly change details like the value of Mc, the main results - i.e. the 
existence of a cutoff and a rather large survival probability of the first gravitationally 
bound objects - should persist. 

Here, we will follow the approach of |2Hj and assume that the microhalos trace the 
mass distribution except for regions very close to the galactic center, where strong tidal 
forces and/or interaction with stars are expected to have destroyed all substructure by 
now. A very simplified model to describe the dumpiness of any halo distribution then 
contains three parameters; the disruption radius Rd within which all substructure is 
completely disrupted, the fraction / of the dark matter outside this radius that comes 
in clumps, and the dimensionless parameter 

1 /d3rei(pci(rd))^ 
po / d3rciPci(rci) 

which describes the effective contrast between an average dark matter clump and the 
local halo density. The numerical results indicate microhalo profiles pd oc r~"' with no 
further substructure and slopes 7 in the range from 1.5 to 2 down to the resolution 
limit of the simulation Taking the average slope to be 1.7, one can in principle 

obtain density contrasts up to 5 ~ 10^ in the most extreme case, depending on how 
far one trusts the extrapolation. (Here, we have assumed the existence of a minimal 
radius inside which the density pd is constant and the inverse of the self-annihilation rate. 
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[{<^a.nn'v) ^(1) Pc\/f^Bw] ^) GQuals the typical time-scale for the formation of the singularity, 
(Gpci)"^, where pci is the average clump density at collapse |39]). 

To calculate the collective contribution from clumps to the total annihilation flux, 
one just has to replace 

(pcDM(r))^ fSpo pcDM(r)6(r - Rd) (8) 

in the source function The only assumption that enters here is that the number 
density of clumps in the galaxy, tracing the halo profile, is essentially constant over 
regions with the size of an individual clump - which should be a good approximation 
for almost any scenario. The simple recipe (jH)) therefore applies not only to antiprotons, 
but to any final state of dark matter annihilations. 

With a clumpy halo distribution, any contribution of dark matter annihilation 
products to cosmic rays will be enhanced. Therefore one must be cautious not to 
violate observational bounds on the flux; this puts severe constraints on the allowed 
parameter space of the clumpy halo model described above. For particles with a long 
free pathlength, however, the by far greatest flux is still expected from the very dense 
region around the galactic center where there is no enhancement due to the presence of 
clumps. To illustrate this, let us consider the case of photons, where the flux observed 
with a detector of angular acceptance AQ is proportional to a line of sight integral 
towards the galactic center, J = J^^dfl Ji^^dl pcDM(r)^- The contribution Jd from the 
clumps alone is therefore given by 

where J = '^Ipcdm^pcdm Pcdm is the smooth halo profile parametrized by ((Zj). An 
NFW profile with = 1 kpc, for example, gives bj^ ^ 6 ■ 10~^, which means that 
the continuous gamma ray signal from KK dark matter annihilations ^2] in that case 
does not overpredict the data taken by HESS jlHI (or EGRET jH]) for clump density 
contrasts as high as 5 < 10 V/ {6 < 10 V/). 

5. Propagation through the galaxy 

For the propagation of antiprotons through the galaxy we use the semiempirical two- 
zone cylindrical diffusion model of ref. JH] to determine the expected number density 
Np[T,r) of antiprotons from LKP annihilations, given the source function In this 
approach, the galaxy is split into two parts and described as a thin disk of radius 
Rh = 20 kpc and height 2hg = 0.2 kpc that contains the interstellar gas (with number 
density Ug = 1 cm~^), and a halo of the same radius and height 2hh (which is not equally 
well known). To simplify things, we assume the same diffusion coefficient 

D{n) = Do (i + ^) (10) 

for both parts, with a dependence on the particle rigidity TZ (i.e. the momentum of 
the particle per unit charge) as indicated. Antiproton losses due to inelastic scattering 
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with the interstellar gas, primarily hydrogen, are taken into account by a corresponding 
loss term in the diffusion equation. Given the cylindrical symmetry, and the boundary 
condition that antiprotons can escape freely at the border of the propagation region, it 
is then possible to find an analytical solution to the problem; its explicit form, expressed 
as an expansion in Bessel functions, is given in pT)] . 

One may refine this simple model by allowing for the presence of a galactic wind, 
blowing particles away from the disk, or by taking into account the effect of the (time- 
varying) solar activity on charged particles near earth. Another possible effect is 
reacceleration by stochastic magnetic fields during propagation or shock waves from 
supernova remnants. In the very high-energy regime that we are interested in, however, 
these effects are of minor importance (see [THl [T7] and references therein) and we do not 
consider them here. 

The allowed range of propagation parameters can be constrained by analyzing the 
flux of stable nuclei, mainly by fitting the observed boron to carbon ratio B/C. This 
has been done in detail for a very similar propagation model |12I and in the following 
we will take a choice of parameters that results in a median flux of antiprotons j43j : 
hh = A kpc, 5 = 0.7 and Dq = 3.5 ■ 10^'' cm^ s~^. The cutoff scale is set to TZq = 1 GV; 
again, its exact value does not influence the antiproton flux at high energies. 

6. Expected antiproton fluxes 

Having solved the diffusion equation for the number density Np{T, r), the expected flux 
of antiprotons from LKP annihilations is readily obtained as 



In addition to this primary source of antiprotons, there is a background contribution 
to the flux that is dominated by antiprotons from inelastic cosmic ray collisions with 
the interstellar medium, mostly protons colliding with hydrogen and helium. Since the 
spectra of both cosmic ray protons and helium nuclei are quite well determined, one 
can again use the B/C ratio to determine the propagation parameters describing the 
diffusion model and obtain a prediction for the background flux that fits the available 
antiproton data remarkably well |44j . 

This background flux is shown in Fig. |21 together with the primary component 
from LKP annihilations. As one can see, the dependence on the halo profile is much 
smaller than for cosmic ray photons, where the induced uncertainties in the flux span 
several orders of magnitude jl?! |M] ; this can easily be understood by noting that the 
halo distributions in question mainly differ very near the galactic center - whereas 
most antiprotons simply do not reach that far through the diffusive halo § . A much 

§ Note, however, that we have adopted here a rather conservative version of the Moore profile and 
assumed a constant density inside a core radius of 0.1 kpc, which is roughly the resolution of the 
underlying N-body simulations. Extrapolating the r~^'^ behaviour down to 10"^ kpc, where higher 
densities are prevented by LKP self-annihilation |39j . gives an enhancement of 5.3 (for low energies) up 
to 10.5 (for high energies) as compared to the NFW case . 
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Figure 2. The antiproton background flux as determined in extrapolated for the 
energy range from 100 GeV to 1 TeV, is shown together with the data taken by the 
balloon-borne experiments BESS 03 and CAPRICE The dashed lines show the 
contribution from LKP annihilations, with an LKP mass of m^ci) = 800 GeV and, 
from bottom to top, an isothermal, NFW and Moore halo profile, respectively. 



greater astrophysical uncertainty is connected to the fact that the B/C analysis actually 
allows a rather large range of propagation parameters; even though this results in 
nearly degenerate background predictions, the flux from primary sources can vary 
almost an order of magnitude at higher energies (for energies below about 1 GeV the 
induced uncertainy spans two orders of magnitude) |13l Ej • A thorough analysis of the 
uncertainties related to the choice of propagation model, however, is beyond the scope 
of this paper and therefore we have adopted a choice of parameters that results in a 
median flux of primary antiprotons, as specified in Section Finally, the uncertainty 
in the local halo density po induces another factor of about 4 for the primary flux. For 
completeness. Fig. IHl shows the dependence of the primary flux on the LKP mass. 

Comparing the expected background and primary fluxes, prospects for the indirect 
detection of KK dark matter through its annihilation into antiprotons clearly do not 
look very promising. The situation changes drastically, however, if one allows for clumpy 
halo distributions. Making use of (jH)), the contribution from the clumps alone may be 
expressed as 

^i = f5Ki%. (12) 

where $p is the antiproton flux for the smooth halo distribution. The formfactor 
decreases with energy and varies from 0.99 (0.61, 0.35) at T = 0.1 GeV to 0.89 (0.47, 
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0.20) at T = 30 GeV for an isothermal (NFW, Moore) profile; above this, it stays 
approximately constant. As an example, Fig. E] shows the resulting flux for an NFW 
profile with f5 = 200 and an LKP mass m^ii) = 800 GeV (a Moore profile with f6 =190 
would give an almost identical spectrum). For such boost factors, an interesting spectral 
feature thus appears at energies right above the currently available data, that can clearly 
be distinguished from the background. As also indicated in the figure, the soon to be 
launched satellite experiment PAMELA [IHl, with its energy range up to 190 GeV, as 
well as the AMS-02 experiment jlH], planned to be installed on the international space 
station by 2007, will have enough resolution to see - or to rule out - such a distortion in 
the spectrum. Already now, one may use the existing measurements in the low-energy 
range to put an upper bound on the allowed amount of halo dumpiness for KK dark 
matter that is much tighter than the bound obtained from gamma rays. In the case of 
an LKP mass of 800 GeV, one has to roughly satisfy f6 < 10'^ in order not to overpredict 
the antiproton flux that is measured at energies below 1 GeV - though a careful analysis 
of this bound should of course take into account the much higher uncertainties connected 
to the expected primary flux at low energies. 

As a further example of how the distortion of the antiproton spectrum for high 
energies might look like. Fig. El shows the case of a hypothetical dark matter candidate 
with mass 7 TeV and similar couplings like the LKP, that was recently proposed as a 
possible explanation for the TeV gamma ray signal from the galactic center jT2]- Due 
to the high mass, however, one needs even larger boost factors to see an effect at all. 
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Figure 4. The solid line shows the expected antiproton spectrum for the case of a 
clumpy NFW profile with f5 = 200 and mg(i) = 800 GeV; the dotted and dashed lines 
give, respectively, the background flux and the contribution from LKP annihilations 
alone. The data points are the same as those of Fig. El in addition, the detectional 
prospects of PAMELA 0S| and AMS-02 02] are indicated by displaying their projected 
data after three years of operation (only statistical errors are included; error bars 
smaller than the symbol size are not shown). For AMS-02, only energies above 100 
GeV are considered. 



7. Conclusions 

As a massive vector particle, the LKP offers an interesting alternative to the most often 
studied dark matter candidate, the supersymmetric neutrahno, which is a Majorana 
fermion. In this article, we have studied the contribution to the cosmic ray spectrum 
from LKP annihilations in our Galaxy, focusing on high-energetic antiprotons. 

We found that the expected antiproton signal is dominated by the background 
flux for standard assumptions about (smooth) halo profiles. However, if one allows 
for clumpy halo distributions, which is supported by recent theoretical as well as 
numerical analyses [HHl IHEj, one can obtain a signal at high energies that clearly can 
be distinguished from the background flux already by the shortly upcoming experiment 
PAMELA jlH] ; the AMS-02 experiment ^ will eventually be able to measure the full 
energy range of interest. One might worry that the required boost factors to obtain such 
a signal are uncomfortably high - but it is important to note that existing numerical 
simulations actually do allow for very high values of the density contrast inside the 
clumps (keeping in mind, though, that they do not yet reach the resolutions that would 
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Figure 5. The same as in Fig. 0] but now for f5 =1000 and a hypothetical dark 
matter candidate with mass 7 TeV that was recently proposed to explain the observed 
gamma ray signal from the galactic center |12j . 



be needed to finally settle this question; eventually, one also has to redo the analysis of 
[^Eni for the specific case of KK dark matter As we have pointed out, even very 

high density contrasts need not violate bounds from gamma ray measurements, which 
might be the greatest concern; in fact, much more stringent constraints come from the 
available data on the (low energy) antiproton flux itself. We have also checked that the 
expected signal from positrons ^1], for the masses and boost factors considered here, 
is not in conflict with the available data. Finally, the boost factors due to clumps that 
are required in the supersymmetric case are often even higher [HH] . 

A full analysis of the astrophysical uncertainties, most notably the impact of differ- 
ent propagation models on the results presented here, is beyond the scope of this article 
and left open for future studies. Another aspect well worth studying is the flux from 
nearby individual clumps; this would add to the average effect of small clumps that 
we have considered in this article and thus further enhance the signal. Furthermore, 
it should be of great interest to carefully compare the predicted spectra for the LKP 
and the neutralino, respectively, and see whether PAMELA or AMS-02 would be able 
to distinguish these two cases. In this context it is an important observation that the 
boost factors required for an LKP-induced distortion of the antiproton spectrum au- 
tomatically result in a distinguished peak signal in the positron spectrum even for the 
relatively high LKP masses considered here ^1] - while such a signal would be absent 
in the supersymmetric case due to the helicity-suppressed annihilation of neutralinos 
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into lepton pairs. As an additional remark, the simultaneous appearance of the positron 
peak in the case of a bump-like distortion of the antiproton spectrum would also pro- 
vide a rather firm evidence for the dark matter nature of the latter and disfavour any 
attempt to merely attribute it to a lack of understanding of the antiproton background 
flux. In fact, the positron excess reported by the HEAT experiment ^ if confirmed 
by PAMELA and AMS-02 - might already be an indication for such a feature in the 
positron spectrum. The high-energy window from about 10 GeV to several 100 GeV 
in the antiproton - as well as positron - spectrum will in any case remain a promising 
place to look for signatures from new physics. 
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